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1) Label root voxel with integer 0; 

2) Construct a queue and line up the root in the queue; 

3) If( There is at least one voxel in the queue ) 

Serve the voxel x on the top of the queue: 
For( each of x*s 26-connected neighbor voxel y) { 
If( y in the volume and has not been labeled yet ) { 
Line up the y in the queue; 

/* label the voxel y */ 
Set disx = 999999 ; 

For( each of y's 26-connected neighbor voxel z ) { 

K( z in the volume and has been already labeled with an integer of ) { 

where if (y,z) is 10, 14 or 17 if the Euclidean distance between y and z is 
l,V2^orV3 , respectively; 

) 

If( dist >dj_) \ 

label y with integer disi ; 
dUsi = dj. ; 

} 

) 

) 

} 

X leaves the queue; 

) 

Else { 

end of calculating the distance map. 

) 



